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Abstract: We study a periodically driven (symmetric as well as 
asymmetric) double- well potential system at finite temperature. We show that mean 
heat loss by the system to the environment (bath) per period of the applied field 
is a good quantifier of stochastic resonance. It is found that the heat fluctuations 
over a single period are always larger than the work fluctuations. The observed 
distributions of work and heat exhibit pronounced asymmetry near resonance. The 
heat losses over a large number of periods satisfies the conventional steady-state 
fluctuation theorem, though different relation exists for this quantity 
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I. INTRODUCTION 



Stochastic Resonance (SR) was discovered barely about two and half decades ago, 
yet it has proved to be very useful in explaining many phenomena in natural sciences [1- 
3] . SR refers to an enhanced response of a nonlinear system to a subthreshold periodic 
input signal in the presence of noise of optimum strength. Here, noise plays a con- 
structive role of pumping power in a particular mode, that is in consonance with the 
applied field, at the cost of the entire spectrum of modes present in it. SR, so defined, 
leaves a lot of liberty as to what is the physical quantity that is to be observed which 
should show a maximum as a function of noise strength [4-23]. In other words, no 
unique quantifier of SR is specified. Also, in order that SR be a bonafide resonance 
the quantifier must show maximum as a function of frequency of the applied field as 
well. For instance, in a double-well system, hysteresis loop area, input energy or work 
done on the system in a period of the driving field and area under the first peak in 
the residence time (in a well) distribution are used to characterize SR as a bonafide 
resonance [4- 17, 19-22]. 

In the present work, motivated by recently discovered fluctuation theorems, we 
show that in an overdamped bistable system input energy per period as well as the 
energy absorbed per period by the system from the bath, i.e, the heat, can be used 
as quantifiers to study SR. Also, it is found that the relative variance of both the 
quantities exhibit minimum at resonance; that is, whenever input energy and heat show 
maximum as a function of noise strength (as also frequency), their respective relative 
fluctuations show minimum. This shows that at SR the system response exhibits 
greater degree of coherence. These fluctuations, however, are very large and often 
the physical quantities in question become non-self-averaging. We study some of these 
aspects in the light of the fluctuation theorems in the following sections. The fluctuation 
theorems are of fundamental importance to nonequilibrium statistical mechanics [24- 
46]. The fluctuation theorems describe rigorous relations for properties of distribution 
functions of physical variables such as work, heat, entropy production, etc., for systems 



far from equilibrium regimes where Einstein and Onsagar relations no longer hold. 
These theorems are expected to play an important role in determining thermodynamic 
constraints that can be imposed on the efficient operation of machines at nano scales. 
Some of these theorems have been verified experimentally [4 7-53]. 



II. THE MODEL 

We consider the motion of a particle in a double- well potential V(x) = — ^ + ^ 
under the action of a weak external field h(t) = Asm(ut). The motion is described by 
the overdamped Langevin equation[44] 

dx dU(x) J _, . 

where U(x) = V(x) — h(t)x. The random forces satisfy (£(£)) = and = 
2jk B T5(t -t), where 7 is the coefficient of friction, T is the absolute temperature 
and ks is the Boltzmann constant. In the following we use a dimensionless form of 
equation(l), namely, 

where U(x) = —\ + \ — xh(t), and the external field h(t) = Asm(ut). Now, £(t) 
satisfies (€(t)£(t')) = DS(t — t'), where D = 2k B T. All the parameters are given in 
dimensionless units (in terms of 7, a and b). We consider A <^ 0.25, so that the forcing 
amplitude is much smaller than the barrier height between the two wells. 

Following the stochastic energetic formalism developed by Sekimoto[55], the work 
done by the external drive h(t) on the system or the input energy per period (of time 
tJ) is defined as[21] 

[to+Tu, dh(t) , 

where h(t) is the drive field which completes its period in time r w . The completion 
of one period of h(t), however, does not guarantee the system coming back to the 



same state as the starting one. In other words, x(t + t u ) need not be equal to x(t) or 
U(x, t + tJ) may differ from U(x,t). The work done over a period W p equals change 
in the internal energy AU = U(x,to + tJ) — U(x,to) and heat Q absorbed over a 
period (first law of thermodynamics), i.e, W p = AU P + Q p . Since x(t) is stochastic, 
W p , AU P and Q p are not the same for different cycles(or periods) of h(t). The averages 
are evaluated from a single long trajectory x(t) (eqn(3)). From the same calculations 
one can also obtain the probability distribution P(W) and various moments of W. 
Similarly, appealing to the first law of thermodynamics as stated above we can obtain 
P{Qp) and P(AUp) and their moments, where the subscript p indicates evaluation of 
the physical quantities over one period of the field. Numerical simulation of our model 
was carried out by using Huen's method[56]. To calculate W p and Q p we first evolve 
the system and neglect initial transients. To get better statistics we calculate W p , Q p 
for 10 6 cycles. In some cases we evaluate W, AU and Q over many periods, n, and 
calculate their averages, again, for 10 6 such entities. 



III. RESULTS AND DISCUSSIONS 



The internal energy being a state variable, average change in its value over a period 
AU P is identically equal to zero. Thus, in the time periodic asymptotic state averaged 
work done over the period (W p ) is dissipated in to heat (Q p ) by the system to the 
bath. Thus, (Q p ) can also be identified as hysteresis loop area. As has been reported 
earlier [19-22], (W p ), the input energy per period, shows a maximum as a function of 
D. Fig(l) shows that (W p ) and (Q p ) coincide, thus both the physical quantities show 
SR. Hence, in this case input energy per period, the heat per period or the hysteresis 
loop area can equally well quantify stochastic resonance. However, in this work we 
focus mostly on the fluctuation properties of these quantities. 

The relative variances Rw and Rq of both W p and Q p respectively show minimum 
(fig(2)) as a function of D. It may be noted that even though (W p ) and (Q p ) are 



identical, fluctuations in W p differ from the fluctuations in Q p . The relative variance 
of Q p is always larger than that of W p for all D. It is also noteworthy that the minimum 
value of the relative variance is larger than one. However, the minimum becomes less 
than one if the averages are taken not over a single period of the field but over a 
larger (integral) number, n > 1, of periods. Therefore, in order to obtain meaningful 
averages of these physical quantities in such driven systems one needs to study over 
time scales much larger than one period so that the averages are significantly larger 
than the deviations about them. Also, as n becomes large, the differences between the 
relative variances of W and Q become insignificant (see inset of fig(2)). Importantly, 
in the system under study, this situation (mean > dispersion) can be achieved by 
increasing the duration of averaging time (or the number of periods, n) more easily 
around the value of D where SR occurs. The minimum of relative variance occurs 
just because the mean value is largest there and not because dispersions are smallest. 
However, as the number of periods n is increased the mean value of heat dissipated over 
the n periods (Q n p) ~ n for all n, whereas the dispersion ~ *Jn for large n so that the 
relative variance decreases with n as 4= and one gets a range of D where the averages 
become meaningful. We have observed numerically that Q np behaves as an independent 
variable only when evaluated over a larger number of cycles n as compared to in case of 
W np . For our present parameters approximately Q np is uncorrelated beyond 10 periods, 
whereas W np is uncorrelated beyond 5 periods. 

In fig(3), we have plotted average heat dissipated (Q p )(— (W p )) over a single period 
as a function of frequency. The values of physical parameters are given in the figure 
caption. The figure shows maximum as shown in earlier literature[21]. Thus (Q p ) acts 
as a quantifier of bonafide stochastic resonance. In the inset we give the corresponding 
relative variance of heat and work as a function of frequency. We observe that heat 
fluctuations are larger than work fluctuations at all frequencies. Near the resonance the 
relative variance shows a minimum. It may be noted that minimum relative variance 
of both quantities W p and Q p are larger than one(fig(2) and fig(3)). 



In fig(4), we plot the probability distribution of W p and Q p for various values of 
D. For low values of D (e.g., D = 0.02) P(W P ) is Gaussian whereas P(Q P ) has a long 
exponential tail as in case of a system driven in a harmonic well and with almost no 
chance of a particle going over to the other well of the double-well potential. As D is 
gradually increased rare passages to the other well becomes a possibility and a very 
small peak appears at a finite positive value of W^or Q p ) (e.g., at D = 0.04). As D is 
increased further, P(W P ) and P{Q P ) become multipeaked and the averages (W p ), (Q p ) 
shifts to their positive values. The distributions become most asymmetric at around 
D = 0.12 (where SR occurs) and the asymmetry reduces again at larger D, fig(4). 
When D becomes large (e.g., D = 0.5) the distribution becomes completely symmetric 
and at such high D values the presence of potential hump becomes ineffective to split 
the distribution into two or more peaks. At very small and very large D values P{W P ) 
is close to Gaussian and so does P{Q P ) but with a slow decaying exponential tail. 
In all the graphs, the distribution of P{Q P ) (P(W P )) extend to negative values of Q p 
{W p ). Finite value for distribution in the negative side is necessary to satisfy certain 
fluctuation theorems. Moreover, P{Q P ) has higher weightage for large negative Q p 
than that of work W p . 

It is worth reemphasizing that W and Q behave as additive (or extrinsic) physical 
quantities with respect to the number of periods n and hence (W np ) and (Q np ) increase 
in proportion to n whereas AU, in this case, is an intrinsic physical quantity and 
— — > as n — > oo. This indicates that the distributions P(W np ) and P(Q n p) both have 
identical characteristics as n — > oo. Therefore, the difference between ( V( Wn p^~^n P ) ^ 

and ( ^P ) vanishes as n — > oo. In the recent literature it is shown that the 
distribution P(W np ) over a large number of periods approaches a Gaussian. Also, if 
one considers W p over a single period by increasing the noise strength, P(W P ) ap- 
proaches Gaussian and satisfies the steady state fluctuation theorem (SSFT). SSFT 
implies [26, 34-36, 44-46, 5 1-53] the probability of physical quantity x to satisfy the rela- 
tion P(x)/P(—x) = exp(/3x), where (5 is the inverse temperature and x may be work, 



heat, etc. In fig(5), the evolution of P(Q np ) is shown as n is increased . As n increases 
the contribution of negative Q to the distribution decreases; besides, the distribution 
gradually becomes closer and closer to Gaussian. There is a contribution to P(Q n p) 
due to change in the internal energy AU which is supposed to dominate at very large 
Q making the distribution exponential in the asymptotic regime[34,35,53]. However, it 
is not possible to detect this exponential tail in our simulations. For large n, P(Q np ) 
approaches Gaussian(inset of fig(5)). The Gaussian fit of the graph almost overlaps 

(Q2 ) — (Q ) 2 

and the calculated ratio, "f, " equals 0.99 for n = 25. This ratio is closer to 

-0 (Qnp) 

one, a requirement for SSFT to hold where P(Q) is Gaussian[22, 44,45]. Fig(6) shows 
the plot of ln( p^Q P ^ ) as a function of (3Q np for various values of n. One can readily 
see that slope of ln( p^g p ^ ) approaches 1 for Q <C (Q np ) for large n. This is a state- 
ment of conventional steady state fluctuation theorem. As the number of periods n, 
over which Q np is calculated, is increased, the conventional SSFT is satisfied for Q np 
less than (Q np ) (e.g., for n = 25, SSFT is valid for Q np less than 0.4, for D = 0.16). 
There exists an alternative relation for heat fluctuation, namely, the extended heat 
fluctuation theorem[34,35]. Here, the distribution function obeys a different symmetry 
property for Q ^> (Q np ) for finite n. As n — > oo, (Q np ) — > oo in this limit, and hence 
conventional SSFT holds which has been clarified earlier in linear systems[53]. 

It is further interesting to investigate effects associated with SR in an asymmetric 
double-well potential involving two hopping time scales instead of one as in the symmet- 

2 4 

ric case. We therefore, consider a scaled asymmetric potential V(x) = :z ^- + \ — cx 
driven by the external field h(t). Fig(7) shows the average input energy (W p ) and 
average heat (Q p ) over a single period as a function of D for various values of the 
asymmetric parameter c. From this figure we find that the peak becomes broader 
and lower as c is increased. The peak shifts to larger values of noise intensities for 
higher c. In other words, the phenomenon of SR is not as pronounced [2] as in case of 
c = 0(fig(2)). It is because the synchronization between signal and particle hopping 
between the two well becomes weak because for c ^ 0, the mean time of passage for 



well 1 to well 2 is different from the mean time of passage from well 2 to well I. As a 
consequence the relative variances Rw and Rq become larger as compared to in case 
of c = 0(fig(2)) as shown in the inset of fig(7). 

In fig(8(a)) and fig(8(b) we have plotted probability distribution P(W P ) and P{Q P ) 
over a single period for different values of asymmetry parameter c for a fixed value of 
D — 0.12, A — 0.1 and u — 0.1. As asymmetry increases the probability for particle 
to remain in the lowest well enhances. Hence particle performs simple oscillation 
around most stable minima over a longer time before making transitions to the other 
well. Hence Gaussian like peak near W ~ or Q ps increases as c increases. The 
weight of P(W P ) for larger values of work(positive as well as negative ) decreases with 
increase in c. However, for P(Q P ), its magnitude at large positive and negative values 
of Q p increases as we increase asymmetry parameter. This contrasting behavior can 
be attributed to the larger fluctuations of internal energy AU P as one increases c. This 
we have verified separately. Due to this contribution of AU P for Q p: nature of P(W P ) 
and P(Qp) are qualitatively different. In all cases for fixed asymmetry c fluctuation in 
heat are larger than fluctuation in work. 

In fig(9) and (10) evolution for P(W np ) and P(Q np ) respectively are plotted for 
various values of number of periods n. We clearly observe that as n increases both 
the distributions tend to become Gaussian distributions with the fluctuation ratio 
,iL, = 1, between their variance V and mean (M) as required to satisfy SSFT as 
mentioned earlier. To satisfy SSFT for heat we have to take larger number of periods 
as compared for work. Only in the large n limit contribution to heat from internal 
energy becomes negligible. In the insets of fig(9) and fig(10) we have shown a Gaussian 
fit (with fluctuation ratio equal to one), which agrees perfectly well with our numerical 
data. Conclusions regarding validity of SSFT for asymmetric case for larger periods 
remain the same as for the symmetric case. 

In summary, we find that SR shown by a particle moving in a double- 
well(symmetric) potential and driven by a weak periodic field can be characterized 



well by the heat (Q p ) dissipated to the bath or the hysteresis loop area. It can 
equally well be characterized by the relative dispersion of (W p ) and (Q p ). At reso- 
nance relative dispersion shows a minimum as a function of both D and u. We also 
show that minimum relative variance can be made less than one by taking long time 
protocols of the applied field. For long time protocols distribution P(Q np ) satisfies 
conventional SSFT for P(Q np ) at Q np <C (Q np ) for finite n[53]. We have also shown 
that SR gets weakened in the presence of asymmetric potential and as a consequence 
fluctuation in heat and work become larger. SSFT too is satisfied for both work and 
heat, when they are calculated over large number of periods. 
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V. FIGURE CAPTIONS 



Fig.l: The average input energy (W p ) and (Q p ) as a function of D for u = 0.1 and 
A = 0.1. 

Fig. 2: The relative variance Rw and Rq over one period are plotted as a function 
of D. In the inset the relative variance Rw and Rq over 25 periods are presented. The 
other parameters are same as in fig(l). 

Fig. 3: The mean heat energy (Q p ) is plotted as a function of uo for D = 0.15 and 
A = 0.1. In the inset Rw and Rq over one period are presented. 

Fig. 4: The distribution P(Wp) and P{Qp) over a single period for various values 
of D: 0.02(a), 0.04(6), 0.06(c), 0.08(d), 0.10(e), 0.12(/), 0.16(#) and 0.5(/i). 

Fig. 5: The evolution of P(Q np ) over different periods is presented. In the 
inset P(Q np ) over 25 periods is plotted together with its Gaussian fit f(Q). Here 
D = 0.12, A = 0.1 and uo = 0.1. 

Fig. 6: The plot of ln(P(Q np )/ P(—Q np )) with temperature (3Q np for different 
periods. Only the range of Q np is presented for which the curves are nearly linear. The 
parameters are same as that in fig(5) except that here D = 0.16. 

Fig. 7: The plot of (W) with temperature for various values of the asymmetry 
parameter, c. Inset shows Rw and Rq for c = 0.1 . 

Fig.8: The distribution of P(Wp) and P(Qp) over a single period for c = 0.0, 
c = 0.05, c = 0.1 and c = 0.15. Here D = 0.12,4 = 0.1 and uo = 0.1. 



Fig.9: The evolution of P(W np ) over different periods for c = O.f.In the inset 
P(W np ) over 25 periods is plotted together with its Gaussian fit f(W). Other 
parameters are same as fig(8). 

Fig. 10: The evolution of P(Q np ) over different periods for c = 0.1. In the 
inset P(Q np ) over 25 periods is plotted together with its Gaussian fit f(Q). Other 
parameters are same as fig(8). 
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